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' Abstract 

. We study the possible asymptotically flat perturbations of Robinson- Trautman 

spacetimes. We differentiate between algebraically special perturbations and general 
' perturbations. The equations that determine physically realistic spacetimes with 

angular momentum are presented. 
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1 Introduction 



^ I Robinson- Trautman [|T7l spacetimes (RT) have been very useful for estimating the total 



d I gravitational radiation in the head-on black hole collision |jT2[ 0. In reference [0 
we have applied these geometries to the description of the total energy radiated in the 
head-on black hole collision with equal mass; and it was shown that our calculations 
agree remarkably well with the numerical exact calculations of Anninos et.al.[H. The 
case of unequal mass black hole collision, was treated numerically in reference PI; and 



our technique based on the use of the RT geometries showed again an impressive 



agreement with the exact calculations. If one wants to generalize these estimates to the 
case of the black hole collision with orbital angular momentum it is necessary to consider 
spacetimes with total angular momentum. This is the main physical motivation for this 
work; the application of these calculations to specific physical models will be presented 
elsewhere ITSll. 
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The Robinson- Trautman vacuum solutions are algebraically special spacetimes which 
are characterized by the existence of a congruence of diverging null geodesies without 
shear and twist. This implies the existence of a prefered family of null hypersurfaces; 
which provides a set of sections at future null infinity. The angular momentum calculated 
on these sections is found to be zero. (In a Bondi frame the asymptotic NP quantities 

and (To are zero; as a consequence the angular momentum vanishes independently of 
the definition used. See for example reference and references therein.) In section § we 
present a short review of the RT vacuum solutions. 

It is natural then to consider the algebraically special perturbations of RT spacetimes; 
since they contain the RT family in the limit when the twist goes to zero. 

Algebraically special spacetimes are physically very important solutions since they 
contain the RT family, the Kerr family of solutions and so also the Schwarzschild solution. 
In section |^ we present in full detail the algebraically special vacuum solutions in terms 
of a new frame, the so called 'standard rotating null tetrad' [p^j . 

The algebraically special perturbations of RT spacetimes are governed by rather com- 
plicated equations; however, the equations simplify when one analysis the geometry in 
the far future in the limit of the asymptotic retarded time going to infinity. This is due 
to the fact that the RT geometry simplifies considerably in this limit. It will be shown 
that the total angular momentum of these perturbations is constant in the asymptotic 
future; therefore these spacetimes are not suitable for the description of a perturbed black 
hole with a dynamic angular momentum. Furthermore, the general algebraically special 
perturbations of RT spacetimes turn out to be divergent in the asymptotic future. These 
characteristics do not seem to be appropriate for the modeling of a spacetime represent- 
ing the last stage after the collision of two compact objects with non-zero total angular 
momentum. We are therefore forced to consider general perturbations of RT spacetimes. 

We dedicate section ^ to the presentation of general perturbations of RT spacetimes. 
The radial dependence of the tetrad is explicitly solved; and we also present the intrinsic 
evolution equations at future null infinity, which are used in the calculation of gravitational 
radiation. 

The important issue of gauge freedom is discussed in both cases, algebraically special 
spacetimes (section |^) and general perturbation of RT spacetimes (section ||). 
Summarizing comments are included in the last section ^. 
Throughout the work we will use the GHP|^ notation. 

2 Robinson- Trautman spacetimes 

Robinson and Trautman[l] studied the vacuum solutions containing a congruence of di- 
verging null geodesies, with vanishing shear and twist. The line element of these metrics 
can be expressed [4] by: 



bar means complex conjugate and A is the two-dimensional Laplacian for the two-surfaces 



where P = P(m,C,C), H 



ds 



,2 
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u — constant, r = constant with line element 

-^dC dC; (2) 

where we are using complex stereographic coordinates {(, for the sphere. 

It is usually convenient to describe this line element in terms of the line element of 
the unit sphere; this is done by expressing P as the product P = V{u, (, C)Po{C, (), where 
Pq is the value of P for the unit sphere. If i denotes the vector field that generates the 
null congruence, then £ — du, £{r) = 1, = and £{() — 0. In other words this is the 
coordinate system adapted to the geometry. It is convenient to use the parameterization 
u such that the mass parameter M{u) = Mq = constant. 

Alternatively one can express the geometry in terms of a null tetrad; since, given the 
complex null vector basis m", m", n") with the properties: 

9ab t = -g^ m« = 1 (3) 

and all other possible scalar products being zero, the metric can be expressed by 

Qoh ^£anb + Ua £b -maffib- ffla mb- (4) 

Using the coordinate system x^,x^, x^) — (u, r, (C + C); KC ~ 0) defined above, one 
can express the null tetrad as: 

L = (dul (6) 



with i — 0,2,3 and where the components U and are: 

e^o, f = (10) 

r r 

with 

eo = V2Po V, eo = -ieo; (n) 
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U = rUoo + f/o + — , (12) 
r 

where 

V 1 + 

t/oo = ^, Uo = --Ky, U, = -^^, (13) 

where the curvature Ky of the 2- metric appearing in equation (|^), is given by 

i^y = I V - ^ QyV ByV + V\ (14) 

the leading order behavior of the second component of the Weyl tensor is — 
and 

X° = 1, = 0, = 0; (15) 

and where 9y is the edth operator, in the GHP notation, of the sphere with metric (^. 

In reference [1] it was found that the vacuum Einstein equation can be reduced to 
a parabohc equation for a scalar depending on three variables, the so called Robinson- 
Trautman equation; which in our notation has the form 

-3MoV = V'' V -V"^ d^V Q^V; (16) 

where S is the edth operator of the unit sphere. We refer to a line element with V 
satisfying this equation as a Robinson- Trautman solution. On the other hand, if the RT 
equation is not required, then the solution is no longer vacuum and there is only one 
component of the Ricci tensor different from zero, given by[^ 

^22 - ^2 ' 

where the (RT) is to emphasize the fact that in this case we are using the null tetrad 
adapted to the null congruence. We refer to this as a Robinson- Trautman geometry. 

At this point it is interesting to remark that the spin coefficients a, r, k, a', r' and e 
are zero in the Robinson- Trautman geometries. 



3 Twisting algebraically special spacetimes 

3.1 Null tetrad components of the exact solutions 

As it was indicated previously, from a geometrical point of view, the twisting algebraically 
special spacetimes are the natural generalizations of RT spacetimes, since they also have 
a congruence of null geodesies with vanishing shear and contain the RT metrics as the 
limit of those spaces when the twist goes to zero. 

This kind of solutions have been extensively studied in the literature in terms of 
different tetrads and coordinate systems (for an account of these formulations of the 
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problem see and references therein). The tetrad and coordinate system used here 
do not agree with those ones. This new tetrad corresponds to what we call "standard 
rotating null tetrad"; they are defined in reference |jl6|, and they can be applied to any 



asymptotically fiat spacetime. In this section we will express the equations governing 
twisting algebraically special spacetimes in terms of the "standard rotating null tetrad" . 

The null tetrad components defined in equations (|^) to @ can be expressed by the 
three scalars M{u, (, (), V{u, (, () and L{u, (, (), and are given by: 

fO f2 f3 

e = ^^, e = ^^, e = ^^, (18) 



r — ic r — ic r — ic 



where 



(.l^LV, {„2 = y2P„V, il = -iil P„ = iJ-^; (19) 



I/ = rC/.o + C/o+ 'f,!"'' -(^ + ^)- , .. . (20) 

[r + ic) [r — ic) \r — ic r + ic I [r + ic)[r — ic) 



where 



Uo, = ^. Uo = -^{Kv + Kl), (21) 



Kv = 2VmV-2WW + V'^, (22) 



Kl = 2LLV'^ + \^LLV + LLV - 2WL + V^L + V^L - 2W Lj V 

+ V'^LL + 2VVLL - 2V^LL + V^LL 

+ 2W V{L + L) + V^{dL + BL) (23) 

/i =dLV'^ L + 2cLV'^ L + cLLV'^ + cSL 

+ cci + ^LcV^L + LdcV^ + ^LdcL-^V^L + Lc^iL-^ 

+ ^LcLV^ + ^LdcV^ + B{c)LV^ + ^BLcV^ + mcV^ (24) 

+ ^QLQcL-W^ + QcV^L + ^QLcV^ 

- ^QcdL L-^ + cdci L-^ + c{Kv + Kl) 

Too = ~LV, (25) 
To = 2icToo - iV^LC, (26) 
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where the operator acting on any function f — f{uX,0 defined by: 
and 

The spin coefiicients corresponding to this tetrad are: 

-1 

P 



r + ic 

(7 = 0, 

K = 0, 



Too , To 

T = 



r + ic {r + ic){r — ic)^ 



r — ic {r + ic){r — ic) (r — ic)'^{r + ic) (r + icY{r — ic)^ 



r + ic (r + ic)^ (r + ic)^ ' 



K = : U{t + T), 

r + ic ^ ' 
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r — ic r + ic {r + ic){r — ic) 2(r + ic)^ (r + ic)^(r 

/5o 



r — ic 



where 



/3' = /9 + t'; 



P() = ^^^^ - LL + ^zc - \ + g^L) T/^ 



To = -To, 



3 V2Vm-^B) + S^V 
Po — 7, ■ 



The Weyl tensor has the properties: 



*2 = 



(r + ic)^ ' 
where 

The components ^3 and ^4 are too comphcated to be presented here. 



The vacuum equations require several relations among the main functions. The equa- 
tion $12 = imposes the condition: 



dL{M + = -3L{M + ifi). (54) 
The other relation has a compact form if one uses the function z{u, (, () defined by 

l. = | = i; (55) 
using this, the equation <l>22 = imposes the relation 

z 2 du du du 

It is interesting to note that using z the quantity fi acquires the compact form 

iz"^ - - - - 

1^ = — (SlSlSlSl - (57) 

It is easy to see that this tetrad goes to the RT tetrad when L goes to zero. 

The representation of the twisting algebraically special spacetimes in terms of standard 
rotating null tetrads, we have just gave, is different from that appearing in the book of 
exact solutions by Kramer et.al. ^ and its references, because we use a different tetrad 
frame. 

The twist of the null congruence, and also the angular momentum of the spacetime, 
are proportional to the quantity c defined above. 

The complex equation (^i)) and the real equation ( ^61) constitute three equations for 
the real functions z and M and the complex function L; in other words one has three real 
equations for four real scalars. In order to understand this under determined system it is 
important to mention the gauge freedom available. 

3.2 The gauge freedom 

These solutions are asymptotically flat. The retarded time null coordinate u defines, at 
future null infinity, a set of sections by the equations u = constant. The main gauge 
freedom can be associated to a different choice of these sections; therefore let us consider 
a transformation of the form 

u = uo{%C,0 + 0{-), (58) 



+ ro(M,C,C), (59) 



w(m,C,C) 



C = C + 0(-); (60) 
r 
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one could also generalized the last equation to a transformation among the sphere coor- 
dinates {(, (), but they are not important for the following equations and therefore we do 
not consider them main gauge freedom. 

The null character of the coordinate u requires 

w = uq. (61) 

Although the vector i is proportional to i, the other tetrad vectors rotate under this 
transformation. From the study of the rotations one deduces that 

V 

V=-, (62) 
w 



V{L - mo) = VL (63) 

and 

M + ^^^=^. (64) 

Let us observe that there are several possibilities associated to the three equations 
(|6^), (|63|) and (|6^ . If the twist of the null congruence is different from zero then one can 
choose uo{u, (, C) so that V = 1, for example. 

If we write the spin weight 1 quantity L in terms of the spin weight quantity J by 
the relation 

L = idJ = imeJ - mmJ; (65) 

where the complex scalar J can be expressed in term of its real and imaginary parts 3?e J 
and QmJ, so that J = 3?e J + i'^mJ; then equation (p5|) can be written in the form 



idyj — dyUo = iSy J; (66) 

where 9y is the edth operator for the sphere with metric dS"^ = y2p2 dC d^. From this 

equation it can be seen that in principle one could use uq to make QmJ vanish; but it is 
not clear that the needed Uq would satisfy that Uq > 0; as it is required for a non-singular 
transformation. If QmJ where zero then one can see that a transformation of the form 
u = u — 7(C, C) would do the job. It is conceivable that for small time variations one 



can use techniques similar to those of references |]TT[ and to prove local existence of 
transformations that make S^mJ = 0. 

Another possibility is to use the freedom involved in the choice of uq to make M to 
be a constant. This can always be made since from equation (|6^ one always solve for mq, 
since the function M, for physically realistic spacetimes, is assumed to be positive definite. 
This choice shows that from the four original unspecified real scalars, one of them is pure 
gauge; therefore the complex equation (pif ) and the real equation ( ^61) constitute a system 
of three equations for three physical real scalars. 
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3.3 Algebraically special perturbations of RT spacetimes 

Having just presented the exact algebraically special solutions in terms of a new frame, we 
proceed here with our program by calculating the equations governing the algebraically 
special perturbations of RT spacetimes. We will also specialize the discussion to the 
asymptotic future regime for the retarded time going to infinity, where the equations are 
solved. 

When one studies the case of a hnear perturbation of the RT hne element, one con- 
siders: 

V^Vrt + Vl (67) 

and 

M = Mo + Ml (68) 

where Vrt satisfies the RT equation, Mq is the constant RT mass and Vl and Ml rep- 
resents the departure from the RT metric that take into account the presence of angular 
momentum; which are zero when L vanishes. 

Prom equation $12 = it is obtained, in first order 

= - 3 LMo - 5^1 8 Vrt - 3 5^1 5 Vrt § Vrt V^^ - 3 Vrt 

- 5Vi?T L V^^ + 9Vi?T S=^L \/^^ - 3 5 L 8 Vrt V^Jt 

- 6 gZ Wrt ^Wrt \/kt - 2 S I BgVflT V^^ - S M^,, - 3 Vrt SSL V^^ 

+ 3gVRT v^RT - 3 Wrt gg^i v^^ + 3 Wrt ssl Wrt v^^ 

+ 6 Wrt ^Wrt BL V^j, + 3 Wrt B'^dL V^j. - 3 Wrt BL V^t (gg) 

- 3 Wrt Wrt L V^t - ^SS'^ V^t + §9'^ 9 Vrt ^Jt 
+ 2 BS Vrt B L + 3 BSL mVRT V^T - \bql V^t 

- 359VRr L V^t + ^§'9'^^ K^r 
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while from the equation $22 = the following relation is obtained 
= - 2 L Wrt B VijT + 2 B'^Wrt Vrt + ^ S^iJT V^t - 2 ^ 9 Vrt BVrt 

+ 2 L ggV^T V^^T - Vrt + ?>VlMo- Vrt V^Jt 

- 12 ^ Mo + 3 Vrt + 2 5L BSV^t V^t + -B^ ^kt + BKrt BSZ 

Vrt 2 

+ g^RT B^L V^2^ - 2 gVRT BVht + BVrt V^t + BVjjt Bl 

+ g^L BK^t - BVht ^Jt - 2 g Vht B\/ht Z 

- gVflT BVi v^j^ + g b v^t Vl v^t + b V^t 

+ 2 gi Bg^jjT vL + qvrt Bgi y^^. + gvRT B^l + 2 Bl ggvRT 



(70) 



'i?T ^ ^^KT y RT ^ Kr ^ ^ y rt "t ^ ^^vrt y rt 

+ gvRT ggL v^^ + Bvrt l y^^^ + bbl Bvrt v^t + 2 BgvRT Bl 



+ -gg^L + 2 gg v.. l v^, + -g^gL vj, + 2 g^gr.. l 
+ g^gVi 

In the asymptotic future; for m — oo it is known @ that 

VRT = l + 0{e~o); (71) 



therefore to first order in the asymptotic expansion of the retarded time, the above two 
equations become 

= -3Li Afo - g - Igg^Zi - ^Bg^i + ^B^g'Li; (72) 

and 

= -Ml, + 3 Mo + + ^Bg^Zi + Ig^gL^ + g2gVi,; (73) 

where the subindex 1 means that we are considering first order quantities in the expansion 
as M — >■ oo. 

Using the complex scalar J which satisfies (p5|), equations (0) and ( |73D become 

= Sg {^mj^ iMo + 3g {^ej^ Mo -QMl.i- S{Q'^S'^^eJi) (74) 

and 

= -Ml, + 3 Vl, Mo - B^g^JJm Ji + B^gV^,; (75) 

where it can be seen from equation ( [7^ ) that if we use the gauge freedom to make M^, 
a constant, then dQmJi must vanish which implies that QmJi = 0. It is clear from the 
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discussion of the previous subsection that by means of a further gauge transformation one 
can make QmJi = 0. Therefore in this gauge the above equations become 

= g (^3 Ji Mo - B'SVi) (76) 

where now Ji is assumed to be real, and 

= 3Vl,Mo + B^S^Vl,. (77) 

The constant M^^ now plays no role, and it can be made to vanish by a further 
reparameterization of the retarded time coordinate. Integrating equation (|76|) we can 
write these two equations in a more comparable form 

3jiMo = 8^5Vi (78) 

and 

-3\4,Mo = 82gVi,. (79) 
It is important to know the spectrum of the operator appearing on the right hand side. 



The eigenvalues of the operator 9^3^ are|]TT| iK(^+i)(^+2) . ^Yieie i is any non- negative 



integer. 

Therefore the leading order behavior for the functions are: 

Ji = JrYUCX) + JTY2miCX)e^ + Oie^o) (80) 

and 

Vl, = VEl r2„(C, C)e"^ + 0(e-^); (81) 

where J{", J™ and are constants and we have used gauge transformations to set to 
zero spurious integrations constants. 

It is observed that Vl^ has the same asymptotic behavior as Vrt] while Ji has a con- 
stant angular momentum term and a possible divergent behavior for the higher moments 
in the limit u ^ oo. 

Therefore these spacetimes are not suitable for the description of a perturbed black 
hole with a dynamic angular momentum. This forces us to consider general perturbations 
of RT spacetimes, that we treat in the next section. 



4 General perturbations of RT spacetimes 
4.1 The vacuum equations 

Let us consider now a general perturbation of the RT spacetimes. We will still work with 
a null tetrad satisfying equations d^-^), with = 0, and the leading order of and S,^ 
also given by eq. (0)! but where now the scalar V is given by 

V = VRT + \Vy, (82) 
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where Vrt is the RT scalar discussed in section A the hnearization parameter and Vx 
the hnear perturbation scalar. Let us emphasize that this means that the tetrad vector 
i°- has no twist, unlike the corresponding one in the rotating null tetrad used in the last 
section. In this section we use asymptotic coordinates and standard tetrad as defined in 
p!6| ; however it is important to notice that we do not use a Bondi system. We use A to 



emphasize in some equations the linearized terms; but it should be understood that it has 
the value 1 at any stage. 

It is convenient to start the study from the vacuum Bianchi identities. 

Assume that in the null hypersurface S, defined by u = mq, of the original RT spacetime 
one is given the function \l/o(uo, 0- Then, noting that \I'o and \I'i are zero in the 
Robinson- Trautman geometries; one obtains, from eq. (3.91) in ref. that, on S, 
is given by 



^1 



j3 



(5vrt ^0 dr'; 



(83) 



where = ^i(mo,C,C)- 

Similarly, from eq, (3.90) of ref. [TU| one obtains that \l/2 is given on E by 



with 



(Mo + A(Mi +i/i)) 



^5) 



where Mq is the constant RT mass. Mi = Mi(uo, C, C) P = p(^0) C; C) are real functions 
and Qvrt is the edth operator of the 2-sphere with the RT conformal factor. 

When calculating \&3 from eq. (3.89) of ref. [|l^], one should take into account the 
first order contribution coming from the edth primed operator; due to the fact that \l/2 
has zero order terms. 

In order to proceed, we need to calculate the first order perturbation of the spin 
coefficients p and a. 

Equation (3.50) of [M is 



dr 



which has solution 



(86) 



(87) 



after appropriately chosen the freedom in the origin of the coordinate r. 



Then eq. (3.51) of [|T0| becomes 



da 
dr 



2pa + 
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which has solution 



(89) 



where cxo = 0-0(^0, CO- 

From eq. (3.52) of one obtains 

dr 



2pr + *i; 



with solution 



(90) 



(91) 



where To = To (mo, CO- 

One can now integrate eq. (3.31) of [|l^ to obtain the radial dependence of It is 
deduced that 



^^ = ^0 + 1 r r'aCdr'- 
r r I 



(92) 



where does not depend on r; then, one can write the integral equation 



r" a ^* dr" ) dr'] 



(93) 



which can be solved by iterations. But, since we should consider only linear terms on A, 
one obtains that in first order 



e = ^ + - / <yildr' 



a a^o^a [' I ,'^^^^dr"\ dr'. 



(94) 



The radial dependence of the components X* is obtained from eq. (3.28) of [T^; which 



IS 



dr 



2f e + 2re; 



(95) 



which implies, in first order, that 



2a 



^0 , 1 



+ — / r"^^irfr" 1 dr' + C.C. 



(96) 



where c.c. means complex conjugate and the Xq do not depend on r; but from the 
coordinate condition at future null infinity [p!6[] , one actually has Xq = 0. 
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In order to determine the dependence of U on the coordinate r, one could use eq. 



(3.26) of |]T0|; which requires to calculate e' first. Equation (3.58) of the same reference is 

-2f/5 + 2r/5'-^2; (97) 



m our case 

9e' 



dr 
therefore 

= -2r iP - P') + 2r (/?' -(3)- (^, + ^^^^ 

= -4Tf - (*2 + ^^2); 

where we must neglect the product involving r; since it is a quantity of higher order. One 
concludes that 

-^ = -{^2 + ^2); (99) 
which implies that U is of the form 

/r pr' 
/ {^2 + ^2)dr"dr'; (100) 

where Uqq and Uq do not depend on r. 

The r dependence of the tetrad components can be explicitly given in terms of the 
potential Wo{u,r,(,C)^ defined from 



by the relations 



^ + Hi 





l(9Wo 

+ 


2 SPVo 








laWo 


2 91^0 









(102) 



with 



[/ = rf/00 + f/o + ^ + ^ + Af/3, (104) 
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where 



and 



where 



_ V 

f/i = -\ i^l + , (105) 



xo = 1, 

X' = >^eo ( + S + ^^^^^ - ^S^V.rWo ] + C.C., 

X' = ( + S + ^^^^^ - ^s^VnrWo ] + C.C, 



To 


2^° 












2^? 











(106) 



To = 5y^j.o-o, (107) 



iTy = '^^v^vV - ^^vV ^vV + v. (108) 

In these equations we are exphcitly denoting the first order terms by introducing the A 
dependency. The scalar V is given by equation (|5^, and therefore one can express Ky 
by 



where 



and 



Kv = Kv^^ + XKv„ (109) 



2 - 2 - 2 

^Vrt = 77 — ^VflT^VflT^^T ~ T72~^VRT^^Jr ^Vrt^RT + V^T) (HO) 



2 

RT 



2 - 2 - 2 

Vrt Vrt vrt 

2Vx - Vx 

+ t73-^VrtVrt ^Vrt^rt + Vx Vrt + TT- Kv^T- (111) 

Vrt Vrt 

The previous relations make the following spinor components of the Ricci tensor to 
vanish 

$00 = 0, (112) 
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^01 = 0, 



(113) 



$11 = (114) 

and 

A = 0. (115) 

The other equations, coming from the other spinor components, involve time deriva- 
tives of the functions. They are: 

d^^) _ IVrtO^Wq ^ 1 fVRTdWp ^ 1^^^ d^Wo 
dudr Vrt dr"^ r-^ \ Vrt dr 2 dr'^ 

+ ^ + Bv^r^VurWo + + ^Moao - 2Kv^^Wo^ , (116) 



^0 ^ 3^^o _ Q^^^^M, + If,) + Sv.,, {Ky^,) ao, (117) 



where /i is actually related to ctq by 



= ^ (S'v^Hx^o - Ql^^ao) ; (118) 
the zero order term of $22 = gives 

-6Mo^ ^ QvRT^VarKvaT, (119) 

which is just another way to write the RT equation, and from the first order one obtains 

- 6Mo-^ = gy«,5y^, Ky^ - 6-^ 3M-^ - Mi ) - 2Mi 

^KT VRT \ Vrt 



VRT Vrt vrt 



~ Tr~^VRT^o + 2-^Qy^^VRTQyR^ao - -^dy^^VRx^yR^ao + Qv^^ctq. (120) 
Vrt vj^ Vrt 

One could also consider the leading order asymptotic expansion of Wq by writing 

Wo-§: + w^, (121) 
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where \1/q = "^qIu, (.(); in this way one has 



^0 = ^ + (122) 



where the second term is of order 0(l/r^). 

In this case equation (|116|) sphts into two equations, namely 



vi/O = 3]^v[/0 + g^^^vl/O - SModo (123) 



and 



dudr Vrt dr'^ \ Vrt dr 2 dr"^ 

\ ( - dWy 1 

+ - ( +Sy«^5y«^lVi + + -Mo do - 2Ky^^W^ 



+ - ^^^Mo. (124) 



Let us notice that equations ( |124| ), ( |123[ ), ( |117| ) and ( |12CI| ) are the evolution equations 



for the first order quantities \I^o, ^5 and Mi respectively; while equation (|119|) 

(the RT equation) is the evolution equation for the zero order quantity Vrt- The scalar 
ao is an arbitrary complex function that affects the radiation content of the spacetime. 
To understand these equations it is important to study the gauge freedoms considered in 
the next subsection. 

4.2 The gauge freedom 

The main gauge freedom admitted in our calculation is of the form 

u = Uoi% C, + ^iMlO + O f 1 y (125) 



+ O (r°) , (126) 



( = C + 0{-]; (127) 



one could also admit a further transformation of the coordinates of the sphere {(, () into 
itself, that was also available in the original Robinson- Trautman geometry; but this will 
not change the following result. 

The condition g"^^ = 1 imposes the relation 

w = iiQ. (128) 
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This asymptotic coordinate transformation is associated to a corresponding null tetrad 
transformation; which in the leading orders is given by 



du = uodu + Mq^ dC + Mqc '^C + ^ ( ~ 



Uqi m m + (J { - 



(129) 



n= — + o(-\=— — + o(- 

du \r J uq du \r 



(130) 



m 



^/2PaV 



w 



r d( \r^ J r 
\/2PqV w V w 



uq ou oQ 



■ n + 



V 



m + O 



(131) 



since the metric expressed in terms of the new null tetrad must coincide with the metric 
expressed in terms of the original null tetrad, it is deduced that 



W Uq 



(132) 



therefore 



m = m T— ^ n + C 



(133) 



The null tetrad transformation equations can be used to write the leading order trans- 
formation relations for the spinor dyad associated to the null tetrad^; namely 



0^ 



Uq O 



rUQ 



(134) 



(135) 



which implies that the regular dyad at future null infinity is given by 



-\Iuq\0 —L 

w " \ ruQ 



—. L 



Uq 



(136) 



I = r = — 



(137) 
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We can now easily calculate the component \l/2 of the Weyl tensor, in leading order, with 
respect to the new null tetrad, obtaining 



where in our case we have 



and 



where 



V, 



Mo + X{Mi+ijl) 
1 



(- [Mo + A(Mi + ijj)] - X QvutI ^Vht^Vut) ; 



:i38l 



^3 = Sv^'o - Byp'o = -I^VutKv^t + 0(A), (139) 



^1 = a'o-2 Uoo cx'o - B^^Uoo = By«,^ + 0(A) (140) 

Vrt 



^'o = A I - (To 1 , (141) 



RT 

and 

P'o = l^v. (142) 
Let us consider a first order transformation generated by 

uo = u~\^{u,C,0; (143) 

then one has 



(144) 



(1 - A7)3 

where it is important to note that the transformation law for do is 

Aao = Aao-Ag2^^^7 (145) 
and therefore, since 7 is real, it can be seen from equation ( |118| ) that 

/i = /X. (146) 

Then, one concludes that 

Ml = 3Mo 7 + Ml + Qv^tI ^v^T^VnT^ (147) 

where it is observed that the left hand side is real, and the first two terms of the right 
hand side are also real; therefore for this equation to make sense one should have the last 
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term also real, which at first sight is not obvious. In order to check this let us note that 
for any quantity H of spin weight s one has 



^v^vH -^v^vH = -sKvH; (148) 
then, taking H = Qyl one has s = 1 and 

SySySyT - ^v^vl = -Ky ^vl', (149) 
so that acting with By on the last expression one obtains 

8ySyBySy7 - BySy7 = -ByKy Sy7 - KyBv^Vl- (150) 

Since dy^v and Bydy are real operators, one concludes that dyKy Sy7 is a real quantity 
for any V. 

Coming back to equation (|147|) we observe that given Mi{u, (() and Vrt{u, CX) it is 
possible to chose 7(m, (, () such that Mi acquires any value that we want. 
Alternatively, from equation ( |132| ) we have 

^ = Vrt + A(Va + 7); (151) 

therefore we could also choose 7 in order to give to Vx any desired value. 

Summarizing, equations ( |123| ), (|117|) and (|119|) are the evolution equations for the 
asymptotic quantities \E'q, \E''j' and Vrt respectively. And, equation ( |120D can be under- 
stood as the evolution equation for Vx if Mi is known; alternatively, it can be understood 
as the evolution equation of Mi if Vx is given. As it was mentioned previously, the func- 
tion (To is free data; which contributes to the description of the structure of the sources 
and to outgoing radiation. 



5 Final comments 



In subsection |3.1| we have presented in full detail the algebraically special vacuum solutions 



in a new frame. The gauge freedom of these solutions is discussed in subsection p^ . 

The algebraically special perturbations of RT spacetimes, studied in section |3.3| , are 
shown to contain divergent behavior in the asymptotic future, and therefore they do not 
seem to provide with an appropriate model for the description of a binary collision. This 
result contrasts with the usual treatments of perturbations of black holes. For example 
in the study of gravitational perturbations of the Kerr geometry, one normally assumes 
perturbations which, with respect to the Boyer Lindquist coordinates, have a t and 
dependence given by e^(^*+'^<^)|[; where E is a constant considered "mostly real and 
positive" 0]. However it is interesting to note that equation (7.286) of reference |§] reads 

C = D - UitM 

where D is a positive real constant, M the mass of the spacetime and if particularized to 
the case of algebraically special perturbation; namely \t'o = 0, then the complex constant 
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C must be zero; which indicates that e*(^*+™'^) has an exponential divergent behavior for 
late times, as we have found. 

One could also ask, how is that this behavior does not appear in the quasi-normal 
modes study of Schwarzschild black-hole; and the answer is that the very definition of 
quasi-normal modes excludes the exponentially growing 53m S < cases^. It is proba- 
bly worth while to remark that the quasi-normal modes can not describe the Robinson- 
Trautman geometry as a perturbation of the Schwarzschild metric. 

In relation to the general perturbations of the RT spacetimes, we have shown in 
section ^ that they provide with a suitable family of spacetimes for the discussion of the 
estimate of the total energy radiated in the collision of two black holes with orbital angular 
momentum. 

In a future work we will apply this construction to the calculations of the mentioned 
estimates. 
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